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Abstract 

It is shown that for four-transmitter systems, a family of four-by-four unit-rate complex quasi-orthogonal space-time 
block codes, where each entry equals a symbol variable up to a change of sign and/or complex conjugation, can be 
generated from any two independent codes via elementary operations. The two independent groups of codes in the family 
generally have different properties of diversity, but the codes in each group have the same diversity provided that the 
differential symbol constellation is symmetric. It is also shown that for four-transmitter systems, an eight-by-four unit-rate 
complex linear dispersion space-time block code can be constructed by using Hurwitz-Radon families of matrices of size 
eight such that diversity three is guaranteed even when all symbols are independently selected from any given constellation. 
This code is so far the only known unit-rate linear dispersion code that has diversity no less than three for four transmitters 
under any given constellation. 

© 2008 Elsevier B.V. All rights reserved. 

Keywords: Space-time block codes (STBC); Orthogonal STBC; Quasi-orthogonal STBC; Non-orthogonal STBC; Hurwitz Radon 
families of matrices; Diversity analysis 


1. Introduction 

Design and analysis of space-time block codes (STBC) for multiple transmitting antennas have been an 
active field of research since the work by Alamouti [1] and that by Tarokh et al. [2]. STBC is aimed to exploit 
the channel diversity between multiple transmitters and multiple receivers to improve the rate of reliable data 
transmission and/or the performance of bit error rate. STBC is also useful for cooperative relays in wireless 
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mobile networks [3-6], where STBC can be used effectively as if between multiple transmitters and a single 
receiver. 

STBC is a mapping (applied at the transmitters) between a sequence of input symbols and multiple 
sequences of output symbols. The number N of the output sequences typically corresponds to the number of 
transmitters. The ratio of the length T of the output sequences over the length S of the input sequence is called 
the rate of the STBC (assuming that both the input symbol constellation and the output symbol constellation 
have the same dimension). The output of the STBC mapping can be denoted by a T x N matrix C(v) where 
the 5x1 vector v represents the input symbols. Assume that the channel following the N transmitters is 
frequency flat, and there are M receiving antennas at the end of the channel. Then, the received baseband 
signals at the destination over a time interval of T symbols can be represented by the T x M matrix Y: 


Y = C(x)H + W, 


( 1 ) 


where H is an N x M channel matrix whose entries may be assumed to be i.i.d. complex Gaussian random 
variables (Rayleigh fading), and W is a T x M noise matrix whose entries may also be assumed to be i.i.d. 
complex Gaussian random variables. With a coherent maximum likelihood decoder, the pairwise error rate 
(PER) P(x x) averaged over the channel fading distributions is upper bounded as follows [2]: 





( 2 ) 


where E s and No/2 are, respectively, the symbol energy and the variance of noise per dimension; the signal- 
to-noise ratio (SNR) may be defined as the ratio of E s over No/2; r is the minimal rank of C{x— x) over all 
possible distinct pairs of the symbol sequences; Vj (j=l,...,r) are the non-zero eigenvalues of 
C(x — x) H C(x — x). To reduce the PER of a code, one must increase r and the minimum of \Yj=\ v j- The 
value of r is called the diversity of the code, and the minimum of YYj=\ v j determines a coding gain. Diversity 
and coding gain are among the key measures of a code. 

A detailed review of STBC is available in [7,8]. For convenience, we will also refer to STBC simply as codes. 
The most attractive codes are perhaps the orthogonal codes [1], which allow the maximum likelihood 
(optimal) detection to be performed independently on each of the individual symbols. But the unit-rate 
orthogonal complex codes exist only for two transmitters [9]. For more than two transmitters, there are only 
fractional-rate orthogonal complex codes [9]. Upper bounds on the rate of orthogonal complex codes are 
explored in [10]. There are also quasi-orthogonal codes that allow the maximum likelihood detection to be 
performed independently on pairs of symbols [11] or even independently on each symbol as shown in [12]. But 
the quasi-orthogonal code given in [11] does not have a full diversity. Various improvements of quasi- 
orthogonal codes are further developed in [13-16]. In [16], it is shown that unit-rate quasi-orthogonal codes 
with maximal diversity products can be constructed by using a finite information symbol set on square and 
triangular lattices. There are also codes that are designed to maximize an orthogonality measure [17]. 
Numerous other codes can be found via [18-24] and the references therein. 

The purpose of this paper is not to present a new code competing against existing ones. But rather, we reveal 
a structural insight into a class of linear dispersion codes whose properties are intrinsically governed by the 
Hurwitz-Radon (HR) families of matrices. We first explore the quasi-orthogonal codes of the type shown in 
[11,25-28]. It will be shown that all 4 x 4 quasi-orthogonal codes, where each entry of the code matrix is a 
symbol variable up to a sign change and/or complex conjugation, can be generated from any two independent 
codes by elementary operations. This result is a fundamental unification of all existing (as well as numerous 
previously unrevealed) 4x4 quasi-orthogonal codes in this category. If all symbols are selected independently 
from a common constellation (which will also be referred to as common constellation condition), the quasi- 
orthogonal codes may have diversity two. More precisely, half of the quasi-orthogonal codes always have 
diversity two under the common constellation condition, and the other half may have diversity either two or 
four under the common constellation condition. If the common constellation is an odd-numbered phase-shift- 
keying, half of the quasi-orthogonal codes actually have diversity four. As illustrated later, both the 2x2 
orthogonal codes and the 4x4 quasi-orthogonal codes can be expressed as linear dispersion codes in terms of 
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the HR families of matrices. This observation motivated us to explore 8x4 linear dispersion codes using the 
HR families of matrices of size eight. It will be shown that a class of 8 x 4 linear dispersion codes constructed 
with the HR families of matrices have diversity three even under the common constellation condition. To our 
knowledge, this is the first 8x4 unit-rate linear dispersion code that is guaranteed to have diversity three when 
all symbols are independently selected from any given constellation. This is an unique insight, which is 
unknown from the previous studies of linear dispersion codes [19,22,29]. A high diversity order regardless of 
symbol constellation is practically useful since it could reduce the physical layer complexity associated with 
constellation constraint. 

In Section 2, we review the HR families of matrices. In Section 3, we show that all 4 x 4 quasi-orthogonal 
codes can be constructed by two independent codes, and their properties are discussed. In Section 4, we 
introduce a class of 8 x 4 linear dispersion codes using the HR matrices of size eight, and show that these 
codes have diversity three under any given constellation. The proof of the diversity three property is a lengthy 
part of this paper. We hope that interested readers will find the proof theoretically insightful as it reveals 
detailed structures in the problem. In Section 5, we provide a simulation example to illustrate the performance 
of the non-orthogonal 8x4 HR code. 

1.1. Notations and terminologies 

• Lower case letters are used for scalars. 

• Upper case letters are used for matrices. 

• Underlined lower case letters are used for vectors. 

• In normal script, * denotes an undetermined quantity. As superscript, denotes complex conjugation. 

• In normal script, j denotes V—1. In subscript, j denotes an integer. 

• Kronecker product is denoted by 0 as defined later. 

• = denotes “equal by definition”. 

• 1 1 is an / x / identity matrix. 

• As superscripts, T denotes transpose, and H denotes conjugate transpose. 

• 91 denotes real part, and 3 imaginary part. 

• Unless specified otherwise, by “symbol”, we mean a complex variable. 

• Unless specified otherwise, two vectors of variables are said to be orthogonal only if they are orthogonal for 
all values of the variables. 

• All other notations are defined the first time they are used. 

2. HR matrices 

2.1. General properties of HR matrices 

Radon Theorem [30]: Within the space of L x L integer matrices, there is a family of p matrices 
{Aq,A\, ..., A p - i} satisfying A 0 = II (the L x L identity matrix) and: 

• Property 1(a): A t Aj = I L - 

• Property 1(b): A- L = —Aj (i> 0). 

• Property 2(a): AjAj = -AjA f (i^j), 

where the maximum value p max of p is governed by L as follows. Let L = 2 a b where b is odd, a = 4c + d and 
0 ^ d ^ 3, then /> max = 8c + 2 d . 

A family of matrices defined above is called an HR family of matrices of size L. The following properties of 
the HR matrices follow readily from Properties 1 and 2(a): 

• Property 2(b): AiAj = -AjAj (i^j). 

• Property 3: For any real vector v, v t AJAj v = S(i — j)\v \ 2 where <5(x) is one when x = 0 and zero otherwise, 
and | v | is the norm of v. 
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• Property 4: Given the matrices Aj,j = 1,..., N, within an HR family, we have 

( N \ ( N \ T N 

U X J A JJ (U x J A jJ = U X P■ 

An HR family of matrices of any size can be constructed from the following 2x2 matrices [30]: 



2.2. HR matrices of size two 


For L = 2, an HR family is {I 2 ,R}. The Alamouti code can be expressed in terms of the 2x2 HR 
matrices, e.g., 


C(sus 2 ) = 


s 1 ^2 

~4 S l 
*l(l) 


vi(2) 


+ j h 


^2(1) 


Jl 


,* 2 ( 2 ) 

*i(2)+jx 2 (l) xi(l)+jx 2 (2) 

v —jci( 1) + jx 2 (2) xi( 2 ) - jx 2 (l) 
where s i = xi( 2 ) + j.v 2 (l) and s 2 = 2 Ci(l) + jx 2 (2). 


Xi(\) 

xi (2) 


+ j^ 


* 2 ( 1 ) 

* 2 ( 2 ) 


(4) 


2.3. HR matrices of size four 

For L = 4, an HR family consists of the following matrices: 

Qo = 3 4 , Qi = P ® R, Q 2 = R 0 1 2 , Q 2 = Q 0 R, (5) 

where ® is the Kronecker product, e.g., 

fa n fli 2 \ B _( anB a ^ B \ 
ya 2 i a 22 J ya 2i B a 22 B J 

The HR families of size four are closely related to 4 x 4 quasi-orthogonal codes as discussed later. 

The following theorem provides the complete set of HR families of matrices of size four. If the first matrix in 
each family is fixed to be identity, there are total 2 x 2 3 = 16 HR families of size four. 

Theorem 2.1. Any HR family of matrices of size four has either one of the following two possible forms : 

fii ={ eo ;±&;± fi 2 ;± 03 } and ^2 ={g 0 ;±g 1 ;±g 2 ;±g 3 }, 

where Go = Q 0 , G\ = Qi[Q<S> Q\, G 2 = Q 2 [—I 2 0 Q], and G 3 = QfQ®I 2 ). The Q t matrices were defined 
previously. {Note that in this paper, ± in one place should be treated as independent of ± in another place unless 
specified otherwise .) 

Proof. The proof of this theorem requires an exhaustive but finite search, which is tedious but feasible. In the 
following, we provide an outline of the proof. The goal is to show that under Properties 1(a), 1(b) and 2, only 
Q i and Q 2 can be valid HR families. 

Let us first search for all possible 4x4 HR matrices satisfying Properties 1(a) and 1(b). Under Property 1(a), 
each entry of an HR matrix F is either zero or ± 1 and each row of F has no more than one non-zero entry. 
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Under Property 1(b), i.e., F T = —F where FAG, we can write 


/ B\,\ 2 \ 

y i? 2 ,l Bi,2 J ’ 


where B/j is a 2 x 2 block matrix. It follows that Bj t = —Bij for i = 1,2 and Bj 2 = —B 2 ,\. Therefore, B\ \ and 
i? 2,2 are both equal to either the 2x2 zero matrix or ±R , where B\ \ and /? 2,2 may have the same or opposite 
signs. Namely, we have 


either 


B\,\ = B2,2 = 0, 
B \2 = - 52,1 *0 


or 


f * 1,1 = ±R, 

J B 2 ,2 = ±R, 

( B\,2 — - 62,1 — 0 . 


( 6 ) 


It is easy to verify that if = B 2 ,2 = 0, the only possible choices of B i >2 are iLU, ±P, =LR, ±0. 

By an exhaustive, finite and feasible search, one can further verify that any HR matrix F of size four has to 
be one of the matrices in or Q 2 . In other words, the above two sets Q\ and Q 2 contain all possible 4x4 HR 
matrices satisfying Properties 1(a) and 1(b). 

We next consider Property 2(a). It is straightforward to verify that Q\ and Q 2 each satisfies Property 2(a). 
Namely, each of Qi and Q 2 (with a fixed set of ± signs) is a valid HR family. In order to prove that there is no 
other possible HR family, we need to show that for any i,j with 1 </j =^3, Q t and Gj cannot co-exist in one 
family. Indeed, it is straightforward to verify that QjGj / — GjQ t , i.e., Property 2(a) is not satisfied by Q t and 
Gj for any ij with 1 3. □ 


2.4. HR matrices of size eight 

For L = 8 , an HR family is determined by the following eight matrices [30]: 

hJ 2 ®R®l2J2®P®R,Q®Q®R\P®Q®R,R®P®Q\R®P®P,R®Q®h (7) 

which is easy to verify. 

An HR family of matrices of size eight (or size integer-power-of-two no less than eight) have the following 
Properties 5 and 6 : 

• Property 5: Given distinct i, k, m, n, then AjA m AjA n is unitary and symmetric. And the eigenvalues of 
AjA m AjA n are ±1 of differing signs. 

Proof. It is easy to verify that AjA m AjA n is unitary. To prove the symmetry, we apply Property 2: 
a t a a t a —a t a.a t a. 

A j ^n — 

= -Aj n A n AjAi 
= -A T n A m AjAj 
= ^ n A i AlA j 

= ( AjA m AjA n ) T . (8) 

We now prove the eigenvalue property. The symmetry and orthogonality of AjA m AjA n tell us that each of its 
eigenvalues is +1 or —1. We only need to prove that the signs of these eigenvalues are always mixed. Suppose 
that the eigenvalues of AjA m AjA n are all equal to one or all equal to minus one. Then, the matrix is either / or 
-/, i.e., 

AjA m AjA„ = ±I L 
=*AjA„ = ±AJ„Aj 

=>A„ = d -.A t AlAj, (9) 
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where the signs in the above three equations are consistent with each other. Considering another member Aj in 
the HR family where / is distinct from m, n, ij, we have 

A A T _ A.aY A.aY 

= ^a,aJa,„aT 
= =F A,A T m AjAj 

= ±A,AjA m Aj, (10) 

where the signs in the above equations are consistent with each other. From (9), we have 

~ A i A l = tA/AJ A m Aj. (11) 

Since A n Aj = —AiAj t , (10) and (11) imply —AiAjA m Aj = AiAjA m Aj and hence A/AjA m Aj = 0. This 
contradicts the condition of non-zero eigenvalues. Therefore, all eigenvalues of AjA m AjA n are ±1 of differing 
signs. □ 

• Property 6: If AjA k Aj n A n ± AjA t A^A r = 0 where all matrices are distinct HR matrices in one family, 
we have 

A h A kx A m A ^ ± (- l ) EA J A h A l A n = 0, 

where E is an integer and 

\ A h •> A k\ •> A m\ i 5 A j\, A( { , Aq { , A ri ] 

is a matrix series generated by exchanging E pairs of matrices in the following matrix series: 

[ A i, A k , A m , A n , Aj, A t , A q , A r \. 


Proof. It suffices to prove that by switching any pair of matrices in AjA^A^An ± AjA t A^A r = 0, the ± sign 
changes to the sign. 

First, we consider the two matrices exchanged are from either the first term AjA k Aj n A n or the second term 
AjA t A^A r . Exchanging any two matrices A and B that have d other matrices in between is equivalent to 
the following process: (a) exchanging A with its next matrix repeatedly until A is right behind B and then (b) 
exchanging B with its preceding matrix repeatedly until B is in the original position of A. Step (a) undergoes 
d + 1 exchanges of neighboring matrices, and step (b) undergoes d exchanges of neighboring matrices. By 
Property 2(a), the resulting expression has changed its sign 2d + 1 times, and therefore there is a net sign 
change. 

We now need to prove the Property 6 for the case where a pair of matrices exchanged are “crossed” between 
the two terms. We first observe the following equivalent expressions: 

AjA k A T m A n ±A]A'A]A r = 0, (12) 

<=>• AjAiAjA k AlA„ ± AjA,AjA,AjA r = 0 

A]A k A T m A n AjA,AjA r = 0, (13) 

where we should notice that only the leading matrices At and Aj are actually exchanged when we move from 
(12) to (13). We now consider a matrix A (from the first term) that has dA matrices proceeding it and another 
matrix B (from the second term) that has ds matrices proceeding it. In order to exchange A and B, we can 
do the following steps: (a) move A to the front of the first term and move B to the front of the second term, 
(b) exchange A and B, and (c) move B to where A initially was and move A to where B initially was. We see 
that step (a) undergoes the sign change d A + d B times, step (b) undergoes the sign change once, and step (c) 
undergoes the sign change d A + d B times. Therefore, there is a net sign change. □ 
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Given the HR family of size eight shown in (7), it can be verified that 

AJA 6 AJAq = AJA 2 AJA 3 . (14) 

(At this stage, it is not clear whether (14) holds under a more general condition.) Let i,k,m,n,j,t,q,r be a 
permutation of [0,1,, 7]. Then, together with Properties 6 and 3, (14) implies the following properties: 

• Property 7(a): Aj.AiAj n A n = ±AjA t AjA r where the sign depends on the ordering of the indices. 

• Property 7(b): A^A„ = =pA^AfAjA t AjA r where the sign is consistent with Property 7(a). 

• Property 7(c): For any real vector v, v T (AjA,-AjA t AjA r )v = 0. 

3. Quasi-orthogonal codes 

We now present a complete family of 4 x 4 quasi-orthogonal codes (or code matrices) of Type I. A 4 x 4 
quasi-orthogonal code matrix of Type I is defined as such that each entry in the matrix is an element from 
the symbol set and a pair of columns of the matrix is orthogonal to the other pair 

(and the two columns in each of the above pairs are not necessarily orthogonal to each other). Note that each 
± is an independent plus or minus sign, and each superscript ( denotes independently the presence or absence 
of complex conjugation. The above definition of quasi-orthogonal code of Type I was used in [11]. It is 
obvious that if £( 51 , 52 , 53 , 54 ) is a Type 14x4 quasi-orthogonal matrix of the four symbols 51 , 52 , 53 , 54 , then 
numerous Type I quasi-orthogonal codes can be constructed by the following elementary operations: 

C^l,J2^3^4) = PrS(±$,±$MgM£)Pc, (15) 

where (k\,k 2 ,k 3 ,kf) is a permutation of (1,2,3,4), P T permutes the rows and/or reverses the signs of none or 
some rows, and P c permutes the columns and/or reverses the signs of none or some columns. Note that the 
above statement is obvious because none of the operations P r , P c , ± and * changes the quasi-orthogonality of 
(15). While the above statement is obvious, a number of Type I quasi-orthogonal code matrices have been 
introduced in the literature without mentioning the connections among them. For beginners, each of these 
codes appears to be a new one. Even for experts, it was unknown whether all existing codes of Type I are 
related to each other by (15). In this section, we will show that not all Type I codes are related to each other by 
(15), but, however, there are only two groups of Type I codes. The codes in each of the two groups are related 
to each other by (15), but no code from one group is related to any code from the other group. 

One can also extend the family of quasi-orthogonal codes by allowing left or right multiplication of a 
diagonal matrix to (15), which is a simple extension of the Type I codes. There are also orthogonal or quasi- 
orthogonal codes where the entries of the code matrix are non-linear functions of the symbol vector 
{ 51 , 52 , 53 , 54 } and/or the symbol constellation is constrained [31]. In this paper, we will not consider any quasi- 
orthogonal codes other than the Type I quasi-orthogonal codes. 

Two codes will be called independent of each other if they are not related to each other according to (15), or 
otherwise dependent on each other. It is clear that a complete set of quasi-orthogonal codes can be generated 
by all independent codes via (15). But we will show that the number of independent codes is two. Examples of 
such two independent codes are also given. All existing codes of this type will be explicitly expressed in terms 
of the two independent codes. 

3.1. Independent quasi-orthogonal codes 

We show next that there are only two independent 4x4 Type I quasi-orthogonal codes. But first, we have 
the following property about complex vectors: 

Lemma 3.1. Let s = s r + j5 f be a 4 x 1 complex vector in the four-dimensional complex space C 4 where s r is the 
real part and s t is the imaginary part. Define the second 4x1 complex vector as p = M r s r + ]M i s i where M r and 
Mi are unitary integer matrices. Then , 5 H /? — 0 holds for all s in C 4 if and only if (to be denoted by iff) a pair of 
elements in s is orthogonal to the corresponding pair in pfor all s in C 4 and the other pair of elements in s is 
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orthogonal to the other corresponding pair in p for all s in C 4 . (The pair-wise orthogonality is equivalent to that in 
Alamouti code [1].) 

Proof. Taking the real and imaginary parts of s H /> = 0 separately, we have s^ p = 0 iff t\ + t 2 = 0 and T = 0 
where t\ = sj:M r s r , t 2 = sj MiS t , and T = sjM r s r — sjs r . 

Because of the independence between s r and s h t x + t 2 = 0 iff t\ = 0 and t 2 = 0. Because both s r and s t are 
any real vectors, we have t\ = 0 iff M r = — Mj and t 2 = 0 iff Mi = — Mj. Similarly, we have T = 0 iff 
M r = Mj. 

The above implies that s H p = 0 holds for all s in C 4 iff M t = —M r and M r = —Mj. 

With the above choice of M r and M if s H /? = 0 is equivalent to sf = 0. We next apply that M r = —Mj 
and each row of M r has only one non-zero entry ±1. Without loss of generality, we can assume that the (/q Jq) 
entry of M r is ±1, i.e., M r (z 0 j' 0 ) = ±1. Then M r (j 0 ,io) = —M r (io,j 0 ) where io^j 0 and all other elements in the 
z’oth andy 0 th rows of M r are zero. This property implies that the z'oth andy 0 th elements in s cancel each other in 
the form ±(s*(io)s*(j Q ) — s*(jo)s*(io))in = 0. In other words, the z'oth and y 0 th elements in s are 

orthogonal to the z 0 th andy 0 th elements in p for all s in C 4 . More explicitly, the z 0 th andy 0 th elements in p are 
±?(jo) and Tfik)- 

Following the same reasoning, s^ p = 0 for all s in C 4 iff the other two elements in s are also orthogonal 
(in the same way as described above) to the other two corresponding elements in p for all s in C 4 . □ 


From Lemma 3.1, the next theorem follows (which corrects a result shown in [32]): 


Theorem 3.1. Define a code set Sq of 4 x 4 quasi-orthogonal codes for the symbol vector s = is\,S 2 9 s^,sfj^, 
where each element in a (normally full rank ) code matrix has the form and two of the four columns in each 
code are orthogonal to the other two over all s in C 4 . Then, the following two codes S\(s\ 9 S 2 9 st„s^) and 
S 2 {si 9 s 2t s^sf)\ 


( Si 

~S4 

S 2 

S 3 ^ 

S2 

S3 

-s* 

~S*4 

S3 

s 2 

~S*4 


^S 4 

Si 

S *3 

5 2 j 



( Si 

S4 

s 2 

5 3 ^ 

S 2 = 

Si 

S3 

-s* 


S3 

Si 

~S*4 

S * 


- S 4 

Si 

S 3 

~ S 2 j 


(16) 


span all the codes in S q via (15), and furthermore S\(s\ 9 s 2i s$,sf) and S 2 (s\ 9 S 2 9 s$ 9 sf) are not related to each 
other via (15). 


Proof. We will say that a pair of codes are dependent of each other if they are related via (15), or otherwise 
independent of each other. Our proof is constructive in that we will construct a largest possible set S of 
independent codes. It is important to stress that permutations of rows and/or columns, permutations of symbol 
indices, change of sign to each row and/or column, and sign and/or conjugation changes to each symbol are all 
variations allowed by (15) among dependent codes. Our proof consists of several steps by which the above 
variations are eliminated from a largest possible set of independent codes. These steps will lead to two possibly 
independent codes -Si and S 2 . These two codes -Si and S 2 are then finally verified to be independent. 

Without loss of generality, we can fix the first column of each code matrix in S to be [si,^,^,^] 7 . 
Furthermore, we can choose S in such a way that the first two columns of each code matrix are orthogonal to 
the last two columns. From Lemma 3.1, it follows that among all code matrices in §, there are no more than 
the following two possible forms up to the variations defined by (15): 

— $2 * ^ 
s\ * 

— ^4 ^ 

*3 */ 

where * (not in superscript) denotes a unspecified entry. Note that in each of T i and T 2 , the first two elements 
of the first column are orthogonal to the first two elements of the third column for all s in C 4 , and the last two 


( Si 

* 

^2 

*) 

1 

( si 

* 

Si 

* 

-*1 

* 

O 

II 

Si 

* 

S3 

* 


* 

S3 

* 

\S4 

* 

*3 

*) 


^4 

* 
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elements of first column are orthogonal to the last two elements of the third column for all s in C. One can 
easily verify that no other possibility exists that is independent of T\ and T 2 . 

Similarly, from T\, one can generate no more than the following four possibilities in S up to (15): 


Ti. i = 


Ti, 3 = 


( Sl 

* 

S2 



( si 

* 

S2 

~ S *A 

S2 

* 

~ S *l 

-sl 

, T h 2 = 

S2 

* 

-s*i 

sl 

^3 

* 

~4 

4 

S 3 

* 

-sl 

4 


* 

4 

4 J 


^s 4 

* 

4 

- 4 ) 

( S\ 

* 

S* 2 


1 

( Sl 

* 

S* 2 

~4\ 

S2 

* 

~4 

~4 

, t ia = 

S2 

* 

~4 

4 

S 3 

* 

-sl 

4 

S 3 

* 

-sl 

-4 

^ S4 

* 

o* 

5 3 

S * J 


^S 4 

* 

4 

4 ) 


Furthermore, from 72, one can generate no more than another four possibilities in § up to (15): 


T 2 ,\ = 


^2,3 = 


( Sl 

* 

-S 2 

-4\ 

1 

( Si 

* 

-S 2 


S2 

* 

4 

4 

» T2,2 = 

S2 

* 

4 

-sl 

S3 

* 

-sl 

4 

S3 

* 

-sl 

4 

^ S 4 

* 

0 * 

5 3 

-4 ) 


^s 4 

* 

5 3 

4 ) 

( S\ 

* 

4 

~4\ 

1 

( Sl 

* 

4 

~4\ 

S2 

* 

4 

-4 

, T 2 ,4 = 

S2 

* 

4 


S3 

* 

-sl 

4 

S3 

* 

-sl 

4 

^S 4 

* 

4 

4 ) 


^s 4 

* 

4 

4 ; 


By filling the second column of each of the above matrices (to satisfy the orthogonality condition), it follows 
that there are no more than the following eight possibilities in S up to (15): 


Tia = 


T 13 = 


T 2 ,\ = 


( Sl 

S 4 

S 2 

-sl) 


( Sl 

s 4 

S2 

S3 

-4 

-sl 

, T u = 

S2 

S 3 

S3 

-S 2 


4 

S3 

S 2 

^ S 4 

Sl 

*3 

4 J 


y S 4 

Sl 

( Sl 

S3 

4 

~4) 


( Sl 

S3 

S2 

-S 4 

_o* 

Sl 

-4 

, T i 5 4 = 

S2 

s 4 

S3 

Sl 

-S 4 

4 

S3 

-Sl 

^s 4 

s 2 

*3 

4 J 


^S 4 

-S 2 

( Sl 

s 4 

_o* 

S 2 

-sl) 


( Sl 

-S 4 

S2 

S 3 

4 

4 

, T 2,2 = 

S2 

S3 

S3 

-S 2 

-sl 

4 

S3 

S 2 

^ S 4 

-Sl 

4 

~4 J 


y S 4 

-si 


-4\ 


S 2 y 


* / 
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/ 1 

Si 

S 2 

-S\\ 

l 

(Si 

Si 

~ S 2 

~ S 4\ 

S2 

-S 4 

s\ 

S 1 


Si 

S4 

^1 

Si 

Si 

-Si 

~ S *4 

S 2 

, 7^2,4 = 

Si 

Si 


S 2 

1 S 4 

S 2 

S 1 

S\ ) 


1 S4 

Si 

4 

*1 / 


If we let S i = r u and S 2 = T i )2 , then it can be verified that 


Th, 3 (A) = PnS 2 (-suS2,S4,s 3 )P 2 i, (17) 

TiM = PnSi(s3,s 2 ,S4,si)P 2 2, (18) 

T 2 ,i(s) = PnS\(s 4 , s 2 , S 3 , s\)P 2 i, (19) 

Tl,2(s) = Pi4S 2 (-SuS 2 ,S3,S4)P 2 4, ( 20 ) 

7"2,3(5) = -Pl5^1 (-^l, *^2, ^4, ^ 3 )^ 25 , (21) 

TiAs) = P\6S 2 (si,s 2 ,S4,S3)P 26 , (22) 

where 

7*11 = P-{\ ,l),+(2,2),+(3,4),+(4,3), (23) 

Pll = P+( 1,1),—(2,2),—(3,3),—(4,4), (24) 

Pl 2 = ^(l,4),+(2,2),+(3,1),+(4,3), (25) 

P22 = 7 ? (l,l),+(2,2),-(3,4),-(4,3), (26) 

P\3 = 7 J +(l,4),+(2,2),+(3,3),+(4,l), (27) 

P23 = 7 5 +(l,l),+(2,2),—(3,4),—(4,3), (28) 

Pl4 = ^-(1,1),+(2,2),+(4,3),+(3,4), (29) 

P24 = 7 5 +(l ) l) ) +(2,2),+(3,3),—(4,4), (30) 

Pl5 = 7 , +(l,l),+(2,2),+(4,3),+(3,4), (31) 

1*25 = ^+(1,1) _(2,2),-(3,3),+(4,4), (32) 

P 16 = 7 > +( l, 1),+(2,2),+(4,3),+(3,4), (33) 

P 2 6 = 7 > +(i ) i) )+ (2,2),-(3,3),+(4,4) (34) 


and P±(i l j l ),±(i 2 j 2 ),±(i 3 j i ),±(uj 4 ) denotes a matrix where the entries at O 2 J 2 ), O 3 J 3 ), O 4 J 4 ) are ±1 and all 

other entries are zero. 

Therefore, there are no more than two independent codes in §, and -Si and S 2 are two possible independent 
codes. 

To prove that S\ and S 2 are indeed independent codes, we need to observe a property from (15). If all 
elements of s come from a common symmetric constellation (symmetric in terms of sign change and complex 
conjugation), then it is obvious from (15) that any two dependent codes Ci(s) and C 2 (s•) must satisfy the 
following identity: 

minrankfCifs;)) = minrank(C 2 (|)), (35) 

5+0 !+o 
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where s is a function of s, satisfying 

f=(dh^ ) ,±^ ) ,±sg ) ,± J W) T . 

In the next subsection, we show that S i and S 2 do not satisfy the property (35). In fact, we will consider an 
equivalent situation where s is replaced by A s and the constellation of A s is symmetric with respect to sign 
change and complex conjugation. Since Si and S 2 do not satisfy the necessary condition (35) as required for 
any pair of dependent codes, S 1 and S 2 are independent. □ 


3.2. Diversity of quasi-orthogonal codes 

The diversity of a code matrix C(s) is the minimum rank of C(A s) over all possible A s ^ 0 where A s is the 
difference between two symbol vectors. For the diversity analysis of quasi-orthogonal codes, we will assume 
that the constellation of A s is symmetric with respect to sign and conjugation. Such a condition is common in 
practice. Because of the property (35) among dependent codes, to study the diversity of all the 4 x 4 quasi- 
orthogonal codes, it suffices to consider the diversity of Sfs) and S 2 (s). 

The diversity of Sfs) is the minimum rank of Df A 5 ) over all A s ^0 where A s = (Asi, As 2 , As 3 , As 4 ) T . Here, 

a(s) fii(s) 0 0 

tt — j^i(s) a(s ) 0 0 

" 0 *) j6 ia) (36) 

0 -j^i(s) a(s) 

b 2 (s) 0 0 

a(s ) 0 0 

0 a(s) b 2 (s) ’ ( ' 37 ^ 

0 Z> 2 (s) a(s ) 

where 


a(s) = ^2\s k \ 2 , (38) 

k=\ 

b\{s) = 23(sisJ + sjs 3 ), (39) 

b 2 (s) = 2iR(s*S4 + S 2 S 3 ). (40) 

It is clear from (36) and (37) that the rank of A(As) is either two or four as long as As ^0. To determine 
the conditions under which A(As) has the full rank (i.e., rank four), it is useful to examine its determinant. 
The determinants of A(s) are given by det(A(s)) = («(y) - bi(s)) 2 (a(s) + bi(s )) 2 or equivalently 

det(A(s)) = [Lk -^ 4| 2 + 1^2 - j^| 2 ] 2 [lj^i +^ 4| 2 + k+j^l 2 ] 2 , (41) 

det(Z) 2 (s)) = [|si + s 4 | 2 + |s 2 + S 3 | 2 ] 2 [|si - s 4 | 2 + |s 2 - s 3 | 2 ] 2 . (42) 


We show next that the conditions for a full rank A (A s) is generally different from that for a full rank D 2 (A s) 
even if the constellation of As is symmetric. (Note that this was needed to complete the proof of Theorem 3.1.) 

Expression (41) would make Lemma 2.1 in [33] more complete and would also enrich Lemma 2.2 in the 
same paper. It follows from (41) that if and only if the set {jAs 3 } and the set {±As 4 } have no common elements 
except zero or the set {As 2 } and the set {±jAs 3 } have no common elements except zero, then A(As) has full 
rank as long as As ^0. 

Expression (42) is given and well discussed in Lemma 2.1 in [33]. It follows from (42) that if and only if the 
set {Asi} and the set {±As 4 } have no common elements except zero or the set {As?} and the set {±As 3 } have no 


Di(s)= 6 i(s) S 1 (s) = 


and 


D 2 (s)=S 2 (s) h S 2 (s) = 


0 

0 

a(s) 

bi(s) 

0 

0 


Please cite this article as: Y. Chang, et al., An insight into space-time block codes using Hurwitz-Radon families of matrices, Signal 
Process. (2008), doi: 10.1016/j.sigpro.2008.02.011 










ARTICLE IN PRESS 


12 Y. Chang et al. / Signal Processing I 1111 / Ill-Ill 



Fig. 1. The set of A s t (crosses) and the set of jA.s' ; - (circles) for 3-PSK. Note that the constellation of A Si is symmetric in terms of sign change 
and complex conjugation, which meets the condition for (35). 

common elements except zero, then D 2 (A s) has full rank as long as A s ^ 0. This is an observation also made in 
[13-16]. 

For example, if the constellation is odd-numbered phase-shift-keying as illustrated in Fig. 1, then the set 
{jAsi} and the set {±As 4 } have no common elements except zero and the set {As 2 } and the set {±jAs 3 } have no 
common elements except zero, and hence Z>i(A 5) has full rank under As ^0. But for the same constellation, 
Z> 2 (As) does not have full rank under As ^0. 


3.3. Quasi-orthogonal codes in terms of HR matrices 

We now illustrate that all quasi-orthogonal codes are linear dispersion codes by expressing the two 
independent codes S\ and S 2 in terms of the FIR matrices. Let the real and imaginary parts of each symbol Sk 
be expressed as Sk = rjc + j h- It is not difficult to verify the following results. For the first code, 



(n 

-r 4 

-r 2 

-r 3 \ 


(h 

-U 

h 

h \ 


7*2 

r 3 

r\ 

-n 


h 

h 

-h 

U 

Si = 

r 3 

~r 2 

r 4 

r\ 

+j 

h 

-h 

-u 

~i\ 


u 

r\ 

-r 3 

r 2 ) 


\U 

i\ 

h 

-h) 


— [0ofi’ — Q\L\, Q 2 l 1 ? — Q 2 H 1 ] + iiQih’ —Chiu —Q\h-> Qoh]■> 

where r_ x = [r x ,..., r 4 ] T and i x = [z 3 , z 4 , — i x , —z* 2 ] T - For the second code, 



( ri 

r 4 

r 2 

-r 3 \ 


(h 

U 

—i 2 

h \ 


r 2 

r 3 

-r\ 

r 4 


h 

h 

h 

-u 

S 2 = 

r 3 

r 2 

-n 

n 

+j 

h 

h 

U 

-h 


\ r 4 

r\ 

r 3 

-r 2 ) 



U 

-h 

h ) 


— [— Q 3 KZ 2 , Q 2 L 2 > QoKr_ 2 , QJl 2 ] + }[—QoKi 2 , Q\i 2 , Q 3 Ki 2 , Q 2 i 2 ]-> (43) 

where r 2 — [—■ r 2 , —n,r 4 ,r 3 ] T , i 2 = [i\, —i 2 , h, — z 4 ] T , and K = —J 2 <8> Q. 
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Table 1 

Expressions of existing quasi-orthogonal codes in terms of two independent quasi-orthogonal codes 


Authors/reference Original code P v Sj P c 


Papadias and Foschini [25] 

hi 

Si S 3 

S4 

\ 

+(h3) 

+( 2 , 1 ) 


(4\ 

+( 1 , 1 ) 

-(2,4) 


4 

“4 4 

~4 


+(4,4) 

+(3,2) 

Si 

S 3 

-(3,2) 

-(4,3) 


S3 

44 4'1 

S 2 




Sl 




^4 

S 3 _ 4 

“4 

/ 




{ 4 ) 



Hou et al. (20) in [26] 

/ 

S 2 S 3 

S4 

\ +( 1 , 1 ) 

+(2,4) 


( Sl 

\ +(U) 

+( 2 , 2 ) 


-Si 

Si -S 4 

S3 

+(3,3) 

+(4,2) 

Si 

-4 

-(3,4) 

+(4,3) 












~4 

4 4 

~ 

^2 



-4 




{-< 

_ 4 4 

4 

/ 



\ s 2 

/ 


Ran et al. (10) in [27] 

(SI 

Si S 3 

■y 4 \ 


+(U) 

+( 2 , 2 ) 


( s '\ 

+(U) 

+(2,4) 


4 

-s\ -s% 

4 


+(3,3) 

+(4,4) 

s 2 

4 

+(3,2) 

-(4,3) 


4 

“4 -4 

4 




4 




\S4 

£3 S 2 

^ 1 / 





\S4/ 



Tirkkonen et al. (10) in [28] 

( Si 

S 2 S 3 

s 4 \ 


+(1,1) 

+(2,2) 


( 

\ +d,l) 

-(2,3) 


—4 

4 ~ 4 

4 


+(3,4) 

+(4,3) 

s 2 

-4 

-(3,2) 

+(4,4) 












S3 

^4 Si 

S 2 





-4 





4 -4 

4) 





\ ‘ S '3 

/ 


Jafarkhani [11] 

( ^ 

s 2 S 3 

s 4 \ 

-(1,1) 

+( 2 , 2 ) 


(Sl 

\ +(U) 

-(2,4) 


—4 

4 — 4 

4 


+(3,3) 

+(4,4) 

s 2 

-4 

+(3,2) 

-(4,3) 


“4 

—4 4 

4 




-4 




\ S 4 

—S3 —S 2 

S\) 




\ 34 

/ 



A column vector of symbols is used instead of a row vector for a more compact form. The arrays of ±(z',y) in the third and fifth columns 
indicate the ±1 entries in the integer matrices P t and P c . Other unspecified entries of P, and P c are zero. 


3.4. Previously published quasi-orthogonal codes 

The previously published 4x4 unit-rate quasi-orthogonal codes of Type I can all be expressed in terms of 
the above two independent codes via (15). Table 1 summarizes these connections. 


4. Non-orthogonal HR codes 

In this section, we present an 8 x 4 unit-rate code and prove that it has diversity no less than three under any 
constellation. Like the 4x4 quasi-orthogonal codes presented in the previous section, this code is also a linear 
dispersion code constructed from HR matrices. This 8x4 unit-rate code is 


C , (a 1 ,v 2 )=MoA 1 + ]A 4 X 2 ,A\X X + jA 5 X 2 ,A 2 X 1 +}A 6 X 2 ,A3Xl + j^7A 2 ], 

^ ^ ^ ^ ^ ^ _ ✓ 


-3 


Va 


(44) 


where x x and x 2 are two real-valued 8x1 symbol vectors, and i = 0,..., 7, are 8x8 matrices from any 
single HR family of size eight satisfying (14). This code is motivated by the structure of a half-rate 4x4 
orthogonal code for four transmitters. In fact, it is easy to verify that the code (44) is an orthogonal code if the 
8x1 complex vector x { + jx 2 is replaced by a 4 x 1 real vector and A,, i = 0,... ,3 are replaced by 4 x 4 HR 
matrices. This code has a very simple and appealing structure. 

The above code, also referred to as non-orthogonal HR code, is a special form of a more general HR code 
introduced in [5], and is also a special form of the linear dispersion codes introduced in [29,34]. 
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One specific form (example) of the non-orthogonal HR code follows from (44) with the HR family 


given in (7): 






/*!(!) 

+ j* 2 (6) 

xi(3)+jx 2 (7) 

xi(4)+jx 2 (8) 

xi(2) +jx 2 (5) \ 


vi(2) 

- j*2(5) 

xi(4) — jx 2 (8) 

—xi(3) + jx 2 (7) 

-xi(l)+jx 2 (6) 


*i(3) 

-j*2(8) 

-xi(l) +jx 2 (5) 

xi(2)+jx 2 (6) 

-xt(4) - jx 2 (7) 


*i(4) 

3 -] x 2 ( 1 ) 

-xi(2) — jx 2 (6) 

-xi(l)+jx 2 (5) 

yi( 3) — jx 2 (8) 


xi (5) 

+ j* 2 (2) 

*i(7)—jx 2 (3) 

xi(8) — jx 2 (4) 

—xi(6) — jx 2 (l) 


xi(6) 

-j*2(l) 

xi(8)+jx 2 (4) 

—xt(7) — jx 2 (3) 

xi(5) — jx 2 (2) 


xi (7) 

- j* 2 (4) 

-xi(5) -jx 2 (l) 

xi(6) — jx 2 (2) 

xi(8) +jx 2 (3) 


^xj(8) 

+ jx 2 (3) 

—xi(6)+jx 2 (2) 

-xi(5)-jx 2 (l) 

-Xi(7)+jx 2 (4) y 


(45) 


where xfk ) is the kt h element of the real symbol vector x /: . 

From the theorem shown next, the code (44) is guaranteed to have diversity three even when all symbols are 
independently selected from any constellation (regardless of the design of the constellation). To our 
knowledge, for four-transmitters systems, the code (44) is the only known unit-rate linear dispersion code that 
is guaranteed to be of diversity (at least) three under any given constellation. This is a useful property in 
practice since any symbol constellation can be applied to this code while the diversity is guaranteed to be no 
less than three. 

Theorem 4.1. Given any x, + jx 2 ^0, the code C(x x ,x 2 ) defined in (44) has a rank no less than three regardless of 
the constellation from which all symbols are independently selected. 

The rest of this section is to prove Theorem 4.1. Since the proof is quite lengthy, we divide the proof into 
several sections as explained next. 

4.1. Outline of the proof of Theorem 4.1 

The proof consists of a sequence of lemmas, and these lemmas are indexed as follows: 

' Lemma 4.1^= Lemma 4.1.1 
Lemma 4.2 <= Lemma 4.1.1 

Theorem 4.11 ( Lemma 4 . 3 .i.i ( 47 ), 

(14). 


Lemma 4.3 <= Lemma 4.3.1 


All lemmas are stated below. The proofs are given in the subsequent subsections in the order shown above. 
Theorem 4.1 results immediately from the three main lemmas: 

Lemma 4.1 (Proof in Section 4.3). The minimum rank of (44) is no larger than three. 

Lemma 4.2 (Proof in Section 4.4). Among any three column vectors in (44), at least two of them are independent. 

Lemma 4.3 (Proof in Section 4.7). Given (14), if any three column vectors in (44) are linearly dependent, they are 
orthogonal to the fourth vector in (44). 

The above three main lemmas are based on the following supporting lemmas: 

Lemma 4.1.1 (Proof in Section 4.2). Given distinct i, j, m, n and the (non-zero) real vectors x x and x 2 , 
the equation AjX x + )AjX 2 = k(A m x x + j A n x 2 ) holds if and only if 

( k=k] +j k 2 = ±j, 

< (^iA^Aj + A n )x 2 = 0, (46) 

I x x = —k 2 Aj A n x 2 , 

where k\ and k 2 are real numbers. There is always a non-zero x 2 that satisfies (AjA^Aj + A n )x 2 = 0. 
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Lemma 4.3.1 {Proof in Section 4.6). Referring to (44), v 2 is orthogonal to v 3 and y 4 if (a) the condition (14) holds 
and (b) there exist two complex scalars k=k x + )k 2 and t=t\ + j ?2 with \t\ 2 + \k\ 2 ^0 and two real vectors x x and 
x 2 such that 

a = fe.£4](^)- (47) 

Lemma 4.3.1.1 {Proof in Section 4.5). The necessary and sufficient condition for (47) to hold is 
I Mx x = 0, 

| X 2 = ~{k\ + l 2 2 )-\A 6 k 2 + A 7 t 2 f(A 0 - A 2 k\ - A 2 ti)x t , (48) 


where 


M = A 2 k 2 + A 2 t 2 + {k 2 + t 2 ) 1 {A 4 — A(,k\ — A 2 t\){A^k 2 + A'jt 2 ) T {Ao — A 2 k\ — A 2 t\). (49) 

One important property of M is 

M t M = c M I + N, (50) 

where cm>0, N is symmetric and orthogonal, and the eigenvalues of N are ±cm of differing signs. The exact 
content of cm is given in (68) and the content of N is given in (68) and (65). Furthermore, M is singular if and 
only if 

[k\+ti=k\+t\+\, 

| t\t2 + k\k2 = 0 . 

4.2. The proof of Lemma 4.1.1 

Given the complex-valued equation A/x x + jAjx 2 = k{A m x x + jA n x 2 ), there are two corresponding real¬ 


valued equations (real and imaginary parts): 

{A i — k\A m )x x + k 2 A n x 2 = 0, (52) 

{Aj — k x A n )x 2 — k 2 A m x i = 0. (53) 

Because A t — k x A m is an invertible matrix for distinct i and m, has to be non-zero. Otherwise, there is no 
non-zero solution for x x or x 2 . From (53), we have 

*1 = k 2 x A T m {Aj - k x A n )x 2 . (54) 

Taking (54) into (52) leads to 

Zx 2 = 0, (55) 

where Z = {A t - k x A m )k 2 l Aj n {Aj - k x A n ) + k 2 A n . It is clear that Zmust be singular, or otherwise there is no 
non-zero solution for ^ 2 - 

With Properties 2 and 4 of the HR matrices, it is easy to verify that 

Z T Z = [k\ + (1 + k\fk 2 2 + 2k\]h + 2AjA m AjA„. (56) 


It is known that Z is singular if and only if Z has at least one zero eigenvalue. Therefore, based on (56) and 
Property 5, the matrix Z is singular if and only if k\ + (1 + k 2 ) 2 k 2 2 + 2k\ = ±2. This equation is equivalent to 
any of the following equations: 

k\ + (1 + k\f + 2k\k\ = ±2kj, 
k\ + 2(k] + \)k\ + (1 + k]) 1 = ±4k\, 

(1 +k\ + k\) 2 = ±4kj. 
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From the above, we see that the minus sign leads to no real solution. After dropping the minus sign, the above 
is equivalent to any of the following: 

(1 . + - 2 k 2 ){\ +*? + *! + 2 k 2 ) = 0, 

(k i + (1 — A: 2 ) 2 )(A:j + (1 + k 2 ) 2 ) = 0 , 

ki = 0 and k 2 = ±l. (57) 

Therefore, k = ±j. Taking this back into (54) and (55) yields the sufficient and necessary conditions on x 2 and 
Xj as shown in the lemma. 

4.3. The proof of Lemma 4.1 

It follows from Lemma 4.1.1 that there are always non-zero x x and x 2 such that any two columns from (44) 
are linearly dependent of each other. For example, if 

Xj = [—1,1, —1,1, —1, —1, — 1, — 1] T , 

x 2 = [-1,1,1,-1,-1,-1,1,1] t 

then the codeword given in (44) has a rank no more than three. 


4.4. The proof of Lemma 4.2 

Based on Lemma 4.1.1 we can now prove that any three columns from (44) have a rank larger than one. 
Suppose that ^4*Xi + jAjX 2 depends on each ot A m x x +)A n x 2 and ^4 ? Xi + jA r x 2 . From Lemma 4.1.1, we have 

{ Xi = —k 2 AjA n x 2 , 

: (58) 

X\ — K 2 /i i A r± 2 

which implies Aj(A n ± A r )x 2 = 0. However, Aj(A n ± A r ) is an orthogonal matrix, which means that x 2 = 0 
and hence Xi = 0. Therefore, AiX x + j Ajx 2 cannot depend on each of A m x x + j^ 4 «X 2 anc ^ A t x x +j^4 r x 2 . 

4.5. The proof of Lemma 4.3.1.1 

Given the complex-valued equation (47), we equivalently have the following two real-valued equations 
(i.e., the real and imaginary parts): 

A 0 x x = A 2 x\k\ + ^ 3 X 1/1 — A(,x 2 k 2 — A 2 x 2 t 2 , 

A 4 X 2 = A 2 x\k 2 + A 2 x\t 2 + A^x 2 k\ + A 2 x 2 t\ 
or equivalently 

04o — A 2 k\ — A 2 t\)x x = —(A 6 k 2 + A 7 t 2 )x 2 , (59) 

( A 4 — A$k x — Ajt x )x 2 = (A 2 k 2 + ^ 43 / 2 ) 2 ^ 1 - (60) 

Recalling Property 4, (A m — A n k x - A/t x ) is always non-singular for distinct m, n and /, and (A m k 2 + A n t 2 ) is 
non-singular for distinct m and n unless k\ + t 2 = 0. So from (59) and (60), k\ + t 2 ^0 unless both Xi and x 2 
are equal to zero. Also from (59) and (60), Xi = 0 if and only if x 2 = 0. 

From (59), we have 

x 2 = -{k\ + t\)~\A 6 k 2 + A 7 t 2 ) T (A 0 - A 2 k x - A 3 t x )x x . (61) 

Also, from (60), we have an equivalent form of x 2 : 

x .2 — (1 + ^1 + t\) ^(A -4 — A(,k x — A 7 t x ) T (A 2 k 2 + A 3 t 2 )x x . (62) 
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Using (61) in (60) yields Mx x = 0 where 

M = M\ + M 2 (63) 

with Mi = A 2 k 2 + A 2 t 2 and 

M 2 = (k 2 + t 2 ) ^(A 4 — Afrki — A 7 t\)(A^k 2 + A'jt 2 ) T (Ao — A 2 k\ — A 2 t\). 

Clearly, x x £ 0 if and only if M is singular. Summarizing the above, we have that (47) holds if and only if M is 
singular, x r satisfies Mx x = 0, and x 2 satisfies (61). 

To reveal a property of M T M, we now apply Property 4 to obtain that M\Mj = (k\ + t 2 )h and 
M 2 M\ = (k\ + 1 (1 + k\ + t\) 2 Is- Furthermore, 

mJm 2 = ( k 2 + t 2 ) l (Ajk 2 + Ajt 2 )(A 4 Ajk 2 + A 4 AJ 12 — k\k 2 I% — A^A^k\t 2 — A 7 Ajt\k 2 — t\t 2 I%) 
x(Ao — A 2 k\ — A 2 t\) 

= M c \ 2 ± + M c 12,2, (64) 

where 

M c 12,1 = ( k\ + t \) \A\k 2 + A^ti) [A 4 AJk 2 + ^4^7 1 2 + A^Aj (— k\t 2 + A:2/i)](yPo — ^2^1 — A^ti), (65) 

Am 

M c 12,2 = - (k\ + tl)~ l (Ajk 2 + Ajt 2 )(t\t 2 + kik 2 )(Ao - A 2 k\ - A 2 t\) 

= - (k\ + tl)~\t\t 2 + kik 2 )(Ajk 2 Ao - kikih ~ AjA 3 t x k 2 + Ajt 2 A 0 - Ajt 2 A 2 k x - t x t 2 h). (66) 
Using Property 2 and (66) yields 

M c \ 2,2 + M J 12)2 = 2(fc[ + 7 2 ) ! (^U2 + k 2 k 2 yi$. (67) 

Therefore, 

M J M = MjMi + M\M 2 + M c12 ,i + M cU , 2 + Mj 121 + Mj 12j2 

— [(^2 + 4) + (^2 + 4) ! (1 + ^1 + t \)“ + 2(A: 2 + t 2 ) l {k\k 2 + CC) 2 ]^ + M c 12,1 + x . (68) 

V_ ^ _/ V_ _ 

cm N 

It is straightforward to verify that 
N = M cm + M T cm 

— 2(kj + <\ f 1 [Ajk 2 A mAo + Ajt 2 A .y/.7o — A '2 ^2 / 1'W , i ^3 — ■ (69) 

\ ^_ ✓ 

W 0 

We note that a m a m = c&I% where = t\ + k\ + (/^C — t 2 k\) 2 . It is also straightforward to verify that 

W 0 W T 0 =c 2 A h. (70) 

We now need to prove Wo^ ± c^/g. Suppose PFo = ±c^/g. Then, 

Ajk 2 A MAo + yf J/221 m^4o — ^2^42^^ ? 1^3 — A 2 t 2 AMk\A 2 = 


which can be rewritten as 

*2^0 - ^3C) = ±f(A M ) T - t 2 Aj(A 0 - k\A 2 ), (71) 

where f(A M ) = AjA 6 k 2 + AjA 7 t 2 + A^A^-k^ + k 2 t\). 

Note that/(d^) is antisymmetric, and/(4 w ) T /(zl^) = c A I&. For A k where k =£(4,6,1), AJA m = f(A M )A].. 
From (71), we have 

{k\ + t 2 2 )A 0 = (A 2 k 2 + A 2 t 2 ){±f{A M ) - AjA 2 (k 2 t , - k it2 )). (72) 
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We now multiply (72) by Aj to yield 

~(k 2 + t 2 )AjA 0 = — (AjA 2 k 2 + t 2 I){±f T (A M ) — AjA 2 (k 2 t x — k\t 2 )} 

= ± [k 2 AyA 2 AjA(, + k 2 t 2 AjA 2 AjA 7 + k 2 (k 2 t x — k x t 2 )AjA 2 A^A 7 ] 

± t 2 ({Am) + {~^)k 2 {k 2 t x — k x t 2 )I^ + AjA 2 (k 2 t x — kit 2 )t 2 . (73) 


Applying Properties 2 and 5, we add Eq. (73) to its transposed version to yield 


\ 


k 2 ±(k 2 AjA 2 AjA6 + t 2 AjA 2 AjA7 + (k 2 ti — k x t 2 )A 3 A 2 AjA 7 ) + (—l)(k 2 t x — k\t 2 )I% 


= 0 . 


/ 


(74) 


Recall kl + tl^O. From the definition of F 0 shown in (74), it is easy to verify that FqFJ = 
[k\ + t\ + (k 2 t\ — t\k 2 ) 2 ]m^{k 2 t\ — t\k 2 ) 2 . Therefore, (74) implies k 2 = 0. Similarly, we can multiply (72) by 
Aj, and then follow the same analysis as shown above to conclude that t 2 = 0. Therefore, we have proven by 
contradiction that Wo 7 ^ ± cjlg. 

Since M J M = cm!8 + N where N has the eigenvalues ±2^ + t 2 )~ l CA, M is singular if and only if 

cm = 2(k\ + t%r l CA ( 75 ) 

which is equivalent to 

\{k 2 + t 2 )~ — (k\ + t 2 + l) - ] - + 4(t x t 2 + k\k 2 ) 2 = 0 
and hence (51). 


4.6. Proof of Lemma 4.3.1 

Given the definitions of v 2 and v 3 shown in (44), it follows that 
v 2 v 3 = (xjAjA^ + xjAjA 6 x 2 ) + (xjAjA 6 x 2 - xjAjA 2 x x )j 

= (xJaJA 6 x 2 - xjAjA 2 x x )j. (76) 

Using (61) and Property 3, the first term (ignoring j) in (76) becomes 
xJaJA 6 x 2 = - (k 2 + tl)~ l x]A]A 6 (A 6 k 2 + A 7 t 2 ) T (A 0 - A 2 k x - A 3 t x )x x 

= — (k\ + t^f)~ x x^ (k 2 A^ Ao — k x k 2 AjA 2 — k 2 t x AjA 3 + t 2 AjA6AjAo 
— t 2 k x AjA 6 AjA 2 — t 2 t x AjA 6 AjA 3 )x x 

= — (k 2 + t^)~ l xl (lAJA^AJAo — t 2 k x AjAoAjA 2 — t 2 t x AjAoAjA 3 )x x . (77) 

Taking (62) into the second term in (76), we can similarly show that 
xJaJA 5 x 2 = (1 + k 2 + t 2 r 1 xjAjA 5 (A 4 - A 6 k x - A 7 t x ) T (A 2 k 2 + A 3 t 2 )x x 

= (1 +k\ + t 2 y l t 2 xJ(AjA 5 A T 4 A 3 - k x AjA 5 AjA 3 - t x AjA 5 AjA 3 )x x . (78) 

Based on (14) and Property 2, we have AJA 6 AJA q = AjA 2 AjA 3 = — AjA 5 AjA 3 . Also, recall 
1 + k 2 + t\ = k\ + t 2 . Therefore, 

vfv 3 = (k\ + t 2 )~ 1 xJ(-/ 2 AjA 6 AjA 0 + t 2 k x A]A 6 A 7 A 2 + t 2 t x A]A 6 A 7 A 3 
— t 2 A^A 3 A^A 3 + k x t 2 A 2 A 5 AjA 3 + tAiAfAsA^ A 3 )x x ) 

= {k\ + t 2 r l t 2 x T l (k x A T l A 6 A T 7 A 2 + t x AjA 6 AjA 3 + k x A T 2 A 5 A T 6 A 3 + I x AJA 5 AJA 3 )x J. (79) 

v_ v _/ 

m n 
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In order to prove 1^3 = 0 , it remains to prove x^Mnx x = 0. By Lemma 4.3.1.1, Mx x = 0 and 
M J M = CMh + N. From the property of N, it follows that range(c ^/8 + N) and range(cA/7s — N) are 
orthogonal complement of each other. Therefore, the solution space of Mx x = 0 is given by the range of 
M 0 =c M h — N, i.e., x x = M 0 v for any real vector v. Then, x\M n x l = 0 if and only if v = 0 for 

any v. The proof of v 7 MqM n Mq v = 0 is straightforward but very lengthy, the details of which are given in 
the Appendix. 

To prove i§y 4 = 0, we need to exchange A 2 with A 3 and Ae with A 2 in the proof for yfy 3 = 0. With 
Property 6 , (14) still holds after the double exchanges. So, the proof of y^v 4 = 0 is basically identical to the 
proof of i§y 3 = 0. 

4.7. Proof of Lemma 4.3 

Each vector in [y x , y 2 ,y 3 ,y 4 ] depends on two out of eight HR matrices satisfying (14). Exchanging any two 
vectors in [vi,v 2 ,v 3 ,V 4 ] is equivalent to exchanging two pairs of HR matrices. With Property 6 , condition (14) 
continues to hold under any even number of exchanges of HR matrices. Therefore, following the same proof 
as for Lemma 4.3.1, if any three vectors in \y x ,y 2 ,y 3 ,y 4 ] are linearly dependent of each other, they must be 
orthogonal to the fourth vector. 


5. Further remarks on the non-orthogonal HR codes 

5.1. Full diversity non-orthogonal HR codes 

Like the quasi-orthogonal codes, the non-orthogonal HR codes can also be made full diversity by 
introducing proper diversity in symbol constellations. While the best method to achieve full diversity of the 
non-orthogonal HR codes is still an open problem, we give one method here to achieve full diversity. Consider 
the codeword as in (45). Let each complex element of x x (i) + j x 2 (i), i = 1,2,3,4, be from the constellation set 
of exp[j 7 i(m /4 +1/ 6 )], m = 1,3,5,7, and each complex element of Xi(0 + ]x 2 (/), i = 5, 6 ,7, 8 , be from another 
constellation set exp(j 7 im/ 4 ), m = 1,3,5,7. Through exhaustive search, it has been verified that the rank of the 
four column code matrix is always four. 


5.2. Simulation 

To illustrate the performance of the non-orthogonal HR codes, we show a simulation example. For a 
system of single receiver and four transmitters, each block of received data can be expressed as 

y = C(x) h + n , (80) 

where we assume 


• The entries of the fading vector h are i.i.d. Gaussian distributed with unit variance. 

• The entries of the noise vector n are i.i.d. Gaussian distributed of variance x 4 x io _01SNR where SNR is 

the dB value of the ratio of the transmitted power over the noise variance. Here, is dependent on 
modulations. For example, <r^ = 2 for 4-QAM (four symbols on the corners of a square of side equal to 2), 

— 1 for QPSK (four symbols uniformly spaced on a circle of unit radius), and o 2 m = 10 for 16-QAM. 
The factor 4 is due to four transmitters. 

• The code matrix C(x) is given by (45). 


An alternative form of (80) is as follows: 


( 


*00 \ 
300 ) 



(81) 
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Fig. 2. Averaged bit error rate versus SNR for four space-time block codes. The data rate for all codes is 2bits/s/Hz. A system of four 
transmitters and one receiver is considered. 


where 


H = 


Elo^m) 

E;=o 


EL>^+4m-)\ 

Em 


(82) 


We used (81) with a sphere decoding algorithm [35-37] to detect x. For each realization of x, we chose an 
independent realization of h and n. 

In our simulation, we compared the four different codes: the quasi-orthogonal code [11], and the full rank 
quasi-orthogonal code with the constellation rotation given in [14,16], the non-orthogonal code (44), and 
the half-rate complex orthogonal code [9]. To ensure the same bit rate, we used QPSK for the first two codes, 
4-QAM for the third code, and 16-QAM for the fourth code. 

The performances of the four different codes are compared in Fig. 2. The non-orthogonal HR code 
(referred to as diversity 3 code in the figure) shows a good performance in a medium range of SNR, i.e., better 
than the original quasi-orthogonal code and even the half-rate orthogonal code. The half-rate orthogonal code 
performs well at very high SNR because of its full diversity. But the full diversity quasi-orthogonal code 
performs the best among the four codes compared. 

It is important to remember that the code (44) has diversity no less than three for any constellation while the 
full diversity quasi-orthogonal code needs to be readjusted for each different constellation. 


6. Conclusion 

In this paper, we have investigated STBC that have strong connections with the HR families of matrices. 
The key contributions are Theorems 3.1 and 4.1. Theorem 3.1 states that the Type I family of all published as 
well as unpublished 4x4 unit-rate quasi-orthogonal codes are simply variations from two independent codes 
shown in (16) via (15). Theorem 4.1 states that the unit-rate code (44) for four transmitters has a rank no less 
than three under any given constellation. To our knowledge, the code (44) is the only known unit-rate linear 
dispersion code that is guaranteed to have diversity (at least) three for four transmitters. This is a useful 
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advantage since it could reduce the physical layer complexity associated with constellation constraint. It 
remains a challenge to discover whether or not there exists a linear dispersion code that guarantees a higher 
diversity than the code (44) for four transmitters over all possible constellations. It is our hope that the 
in-depth analysis shown in this paper will motivate and help this pursuit. 
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Appendix A. Proof of v = 0 

The proof is relatively lengthy. We will repeatedly apply Properties 1-7 of the HR matrices as well as 
condition (14). To help the presentation of the proof, we will use the sum table to be introduced next. 


A. 1. Introduction of the sum table 

For example, to describe the sum ai AjA 2 AjA 4 + oc 2 AJA^AJA 3 , we use the following sum table: 


coefficients 

products of matrices 

Ot\ 

[1234] 

OL2 

[5673] 


where each group of four integers from (0,1,2,3,4,5,6,7) corresponds to a product of four 
matrices, i.e., 

±[ijmn]= ± AjAjAj n A n . (84) 

To describe p l AjA 4 (otiAjA 2 A 2 A 4 + oc 2 AJA6AJA 3 ) as another example, we will use 


coefficients 

Products of matrices 

/5i«i 

Pia 2 

[24] [1234] 

[24] [5673] 


From Properties 1(a) and 2(a), it is easy to verify that [241234] = [13]. Because of AJA^AJAo = 
AjA 2 AjA 3 = AjA 4 AjA 3 , we can write [245673] = [245367] = [1670][67] = [01]. Therefore, we can also 
express 

p 1 aJa 4 (oc 1 aJ a 2 aJ a 4 + oc 2 AjA 6 AjA 3 ) 


as 


coefficients 

Products of matrices and simplification 

(3-iai 

Pia 2 

[24] [1234] = [13] 

[24] [5673] = [245367] = [1670] [67] = [01] 
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A.2. Main body of the proof 

From the definition of Mo, it is straightforward to verify that 
M 0 = cmI 8 — N 

= 2{k 2 2 + ~ 2k 2 {k\ + 1 2 ) 1 Aj[A4Ajk 2 + ^ 4^7 ^2 + A^Aj(—k\t 2 + k 2 t\j\(Ao — A^tx) 

— 2{k\ + 4) _ 1 ? 2^3 + A<\A^t2 + A(,Aj( /^i t 2 + Z^lXK^O ~ k\A 2 ). 

With C ]=^2 + ^ an d C 2 =k\t 2 + ^ 2 ?i, we can write 


( 87 ) 


Mq = 


coefficients 

products of matrices 

2 Cl 

h 

—2&2Cj _1 

[2460] 

-2k 2 t 2 c^ 1 

[2470] 

— 2k 2 c 2 c^^ 

[2670] 

2k%t 1 Ci 1 

[2463] 

2k 2 t 2 t 1 c { 1 

[2473] 

2k 2 tic 2 c f 1 

[2673] 

-2k 2 t 2 c{ x 

[3460] 

-2t 2 t 2 c{ x 

[3470] 

— 2c2^2Ci 1 

[3670] 

2/cifc 2 t2C]“ 1 

[3462] = -[2463] 

2kit 2 t 2 c^ 1 

[3472] = -[2473] 

2k 1 c 2 t 2 c^ 1 

[3672] = -[2673] 


( 88 ) 


Therefore, 


M 0 : 


coefficients 

products of matrices 

2C! 

h 

-2k\cf 1 

[2460] 

-2k 2 t 2 cf l 

[2470] 

-2k 2 c 2 c{ x 

[2670] 

2 ^ 2 ^ic^ 1 — skik^c^ 1 = 2c 2 k 2 c^ 1 

[2463] 

2k 2 t 2 t\c{ x — 2k\t 2 t 2 c[ l = 2c 2 t 2 c[ 1 

[2473] 

2k 2 t\c 2 c { 1 - 2kic 2 t 2 c^ 1 = 2 c 2 C^ x 

[2673] 

—2k2t 2 cf 1 

[3460] 

-2t 2 t 2 cf 1 

[3470] 

-2c 2 f2C]“ 1 

[3670] 


(89) 
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From (79) we have 

Afjv = k\Aj AfiAj A 2 + t\Aj A(,Aj A3 + /C1AJA5AJA3 + t\Aj A5AJ A3 

= fci [ 1672 ] + / i [ 1673 ] + ^ i [ 2563 ] + * i [ 2573 ]. ( 90 ) 

To compute the product MqM n (where M 0 is symmetric), we have 


coefficients 

products of matrices and simplification 

group index in (95) 

2k\C\ 

/ 8 [1672] - [1672] 

1 

^kik^c^ 1 

[2460] [1672] = [4017] 

2 

-2k\k 2 t2c[ 1 

[2470] [1672] = -[4016] 

3 

— 2&i&2C2Cf 1 

[2670] [1672] = -[01] 

4 

2/CiC2 &2C^ 1 

[2463] [1672] = [4317] = [6025] 

because [1670] = [5243] [5620] = [1743] 

5 

2k\C<it2C{ V 

[2473] [1672] = -[4316] = [7025] 

because [1670] - [5243] 44 [5270] - [1643] 

6 

2kic%c^ 1 

[2673] [1672] - -[31] 

7 

-2kikzt2c[ l 

[3460] [1672] = -[3460] [5043] = [65] 

8 

-2kit 2 t2C^ 1 

[3470] [1672] = -[3470] [5043] - -[57] 

9 

—2kiC2t2C^ 1 

[3670] [1672] = -[3012] 

10 


( 91 ) 


coefficients 

products of matrices and simplification 

group index in (95) 

2iici 

Is [1673] = [1673] 

11 

—24&2 c r 1 

[2460] [1673] = -[2460] [5240] = [56] 

8 

- 2tik 2 t 2 Ci 1 

[2470] [1673] = -[2470] [5240] - [75] 

9 

—2tik 2 c 2 c{ 1 

[2670] [1673] m -[2013] » [3012] 

10 

2£ic 2 fc2<h 1 

[2463] [1673] = [2417] = -[6053] 

because [1670] = [5243] 44 [5630] = [1247] 

12 

2tiC2t 2 C 1 1 

[2473] [1673] = -[2416] = -[7053] 

because [1670] = [5243] <S4> [5703] = [1624] 

13 

2t\C2C^ 1 

[2673] [1673] = [21] 

14 

-2tik 2 t2Ci 1 

[3460] [1673] = [4017] 

2 

-2 t 1 t 2 t 2 c^ 1 

[3470] [1673] = -[4016] 

3 

— 2fiC2f2 c i 1 

[3670] [1673] = -[01] 

4 
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coefficients 

products of matrices and simplification 

group index in (95) 

2k\C\ 

7s [2563] = [2563] 

2 

-2/ci/clcr 1 

[2460] [2563] = -[4053] 

1 

- 2k\k 2 t 2 c [ 1 

[2470] [2563] = [2470] [1407] = -[12] 

14 

— 2k\k 2 c 2 c{ 1 

[2670] [2563] = [7053] 

13 

2kic 2 k 2 c^ 1 

[2463] [2563] = [45] 

15 

2k 1 c 2 t 2 c± 1 

[2473] [2563] = -[4756] = -[3012] 

because [1670] = [5243] <^> [1230] = [5647] 

10 

2&iC2C^ 1 

[2673] [2563] = -[75] 

9 

-2kik 2 t 2 c[ l 

[3460] [2563] = -[4025] 

11 

—2kit 2 t 2 c^ 1 

[3470] [2563] = [3470] [1407] = [13] 

7 

— 2k\c 2 t 2 c^^ 

[3670] [2563] = [7025] 

6 


(93) 


coefficients 

products of matrices and simplification 

group index in (95) 

2*ici 

is [2573] = [2573] 

3 

—2fiA:|c]" 1 

[2460] [2573] = -[2460] [1406] = -[21] 

14 

-2t\k 2 t 2 c{ x 

[2470] [2573] = -[4053] 

1 

— 2t\k 2 c 2 c^^ 

[2670] [2573] = -[6053] 

12 

2tic 2 k 2 c^ 1 

[2463] [2573] = -[4657] = [3012] 

10 

2t\c 2 t 2 c{ 1 

[2473] [2573] = [45] 

15 

2t x c\c{ v 

[2673] [2573] = [65] 

8 

—2tik 2 t 2 Ci 1 

[3460] [2573] = -[3460] [1406] - [13] 

7 

— 2t\t 2 t 2 c^^ 

[3470] [2573] = -[4025] 

11 

—2t\c 2 t 2 c^ 1 

[3670] [2573] = -[6025] 

5 


The following table shows how the common terms in the previous four tables are combined. 
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Group index 

Final Sum 

1 

2fcici[1672] 

2 

2/ciCi[2563] 

3 

2t\C\ [2573] 

4 

0 

5 

0 

6 

0 

7 

-2fcic|q 1 [31] + (-2fcit 2 t 2 cr 1 )(-[13]) + (-2t 1 k 2 t 2 cp)[13] 

8 

— 2kik2t2Ci 1 [6b] + — 2£i&2C^ 1 [56] + 2£ic§C:f I 'f65] 

9 

2M2GC, 1 1 5 71 f 2t l k 2 t 2 c^^A -2fcic|cr 1 |75] 

10 

0 

11 

2tiCi [1673] 

12 

0 

13 

0 

14 

2fiC2C^ 1 [21] + 2k 1 k 2 t 2 c{ 1 [l2] + 2G^| c r 1 [21] 

15 

0 


Therefore, 


coefficients 

products of matrices 

2k\Ci 

[1672] 

2/ciCi 

[2563] 

2*ici 

[2573] 

2Gci 

[1673] 

2fciC2C^ 1 + 2kit2Ci X — 2G 

[13] 

2t x c\c{ 1 + 21^1^ — 2 k\k 2 t 2 C^ 1 

[65] 

2k\ c^c^ 1 + 2kit2C± x — 2t\k2t2C± 1 

[57] 

— 2 tiC 2 Ci 1 — 2t 1 h%c^ 1 4- 2 k\k 2 t 2 Ci 1 

[12] 


(96) 


We now need to multiply each term of MqMn, i.e., (96), by Mo from right. Note that 
v T (AjAiAjA t AjA r )v = 0 and v t (AJA;)v = 0 for any vector v and distinct indices. So, we will use [kijtqr]~ 0 
and [ki\~ 0 to denote that the corresponding terms become zero after being multiplied by v from left and right. 
We will use 

E\=2k\c\c\ x + 2k\t\c\ l — 2t\k2t2C^ x , (97) 

E2=2t\c\c\ X + 2t\k\c^ x — 2k x k 2 t2C- x x , (98) 
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2fci ci [1672] x M 0 = 


coefficients 

products of matrices and simplification 

group index 

2kic\2ci 

[1672]I 8 = [1672] 

1 

—2kiC\2k2ci[ 1 

[1672] [2460] = [1740] 

2 

-2k 1 c 1 2k 2 t 2 c^ 1 

[1672] [2470] - -[1640] 

3 

—2kic\2k 2 c 2 c^[ l 

[1672] [26701- —gO] ~ 0 


2kici2c 2 k 2 c± 1 

[1672] [2463] = [1743] 

4 

2kic\2c 2 t 2 c[ 1 

[1672] [2473] = -[1643] 

5 

2kiC\2c2C [f 1 

[1672] [2673] = -[13] ~ 0 


^kic^k^c^ 1 

[1672] [3460] m [172340] ~ 0 


—2k\C\2t 2 t 2 c^ ^ 

[1672] [3470] - -[162340] ~ 0 


— 2kiCi2c 2 t 2 Ci 1 

[1672] [3670] = -[1230] 

6 


(99) 


2^i ci [2563] x M 0 = 


coefficients 

products of matrices and simplification 

group index 

2&iCi2ci 

[2563] I 8 = [2563] 

2 

—2fcici2&2 Ci 1 

[2563] [2460] = -[5340] 

1 

—2k\Ci2k 2 t 2 c[ 1 

[2563] [2470] = -[563470] ~ 0 


—2kic\2k 2 c 2 c^ 1 

[2563] [2670] = [5370] 

7 

2kici2c 2 k 2 ci 1 

[2563] [2463] = [54] ~ 0 


2kici2c 2 t 2 c^ 1 

[2563] [2473] = -[5647] 

6 

2kiC\2c2Ci 1 

[2563] [2673] = [57] ~ 0 


- 2kici2k 2 t 2 c if 1 

[2563] [3460] = -[2540] 

8 

—2kiCi2t 2 t 2 c 1 1 

[2563] [3470] = [256470] ~ 0 


— 2k\ ci 2c 2 f 2 

[2563] [3670] = [2570] 

5 


( 100 ) 
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2tici[2573] x M 0 = 


2£id[1673] x M 0 = 


coefficients 

products of matrices and simplification 

group index 

2£ici2ci 

[2573]J 8 = [2573] 

3 

-2tiCi2£;2d 1 

[2573] [2460] = -[573460] - 0 


— 2t\C\2k 2 t 2 c-± i 

[2573] [2470] = -[5340] 

1 

-2tiCi2A:2C2c] 1 

[2573] [2670] = -[5360] 

9 

2fici2c 2 /c2c) 1 

[2573] [2463] = -[5746] 

6 

2tici2c 2 t 2 c^ 1 

[2573] [2473] = [54] ~ 0 


2fici2c2cj 1 

[2573] [2673] = [56] ~ 0 


- 2t 1 ci2k 2 t 2 Ci 1 

[2573] [3460] - [257460] - 0 


-2tiCi2t 2 t 2 ci 1 

[2573] [3470] = -[2540] 

8 

—2t\C\2c2t2C\ 1 

[2573] [3670] - -[2560] 

4 


coefficients 

products of matrices simplification 

group index 

2£ici2ci 

[1673] Js = [1673] 

8 

—2fici2fc|c^ 1 

[1673] [2460] = [173240] ~ 0 


-2t\C\2k 2 t 2 c{ % 

[1673] [2470] = -[163240] ~ 0 


-2tiCi2k 2 c 2 Ci 1 

[1673] [2670] = -[1320] 

6 

2tici2c 2 k 2 c^ 1 

[1673] [2463] = [1724] 

9 

2t\C\2c 2 t 2 c ^ 1 

[1673] [2473] = -[1624] 

7 

2fici2c2d -1 

[1673] [2673] = [12] ~ 0 


-2tici2k 2 t 2 c^ 1 

[1673] [3460] = [1740] 

2 

— 2tiCi2f2^2 c i 1 

[1673] [3470] = -[1640] 

3 

-2£ici2c 2 t2c]“ 1 

[1673] [3670] = -[10] ~ 0 



( 101 ) 


( 102 ) 
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coefficients 

products of matrices and simplification 

group index 

Ei2 Cl 

[13]/ 8 = [13] ~ 0 


-E 1 2k 2 2 cY 

[13] [2460] = [132460] ~ 0 


-Ei2k 2 t2C^ 1 

[13] [2470] - [132470] ~ 0 


-E x 2k 2 c 2 c{ 1 

[13] [2670] = [132670] ~ 0 


E\2o 2 k 2 Ci 1 

[13] [2463] = -[1246] 

7 

E\2c 2 t 2 Ci 1 

[13] [2473] = -[1247] 

9 

E^c^c^ 1 

[13] [2673] «* -[1267] 

1 

—Ei2k 2 t 2 C^ 1 

[13] [3460] = [1460] 

3 

—E\2t 2 t 2 Ci 1 

[13] [3470] = [1470] 

2 

— E\2c 2 t 2 Ci 1 

[13] [3670] = [1670] 

10 


coefficients 

products of matrices and simplification 

group index 

Ei2c\ 

[57]/„ = [57] ~ 0 


-E x 2klc{ x 

[57] [2460] = [572460] ~ 0 


—Ei2k 2 t 2 Ci 1 

[57] [2470] = [5240] 

8 

— E\2k 2 c 2 c^^ 

[57] [2670] = [5260] 

4 

Erfc^c^ 1 

[57] [2463] = [572463] ~ 0 


Ei2c2t 2 Ci 1 

[ 57 ] [2473] = [5243] 

10 

E^c^c^ 1 

[57] [2673] - [5263] 

2 

—E\2k 2 t 2 Ci 1 

[57] [3460] = [573460] - 0 


—Ei2t 2 t 2 C! 1 

[57] [3470] = [5340] 

1 

—i?i2c2t2Ci 1 

[57] [3670] = [5360] 

9 
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coefficients 

products of matrices and simplification 

group index 

~E 2 2ci 

[12] Jg = [12] ~ 0 


E 2 2k%cY 

[12] [2460] = [1460] 

3 

E 2 2k2t2C 1 1 

[12] [2470] = [1470] 

2 

E 2 2k 2 c 2 E[ l 

[12] [2670] = [1670] 

10 

-Erfc^c^ 1 

[12] [2463] = [1463] 

5 

—E 2 2c 2 t 2 c^ 1 

[12] [2473] = [1473] 

4 

-Ertclc^ 1 

[12] [2673] = [1673] 

8 

Etfk^Ci 1 

[12] [3460] = [123460] - 0 


E 2 2t 2 t 2 c x 1 

[12] [3470] = [123470] ~ 0 


E 2 2c 2 t 2 c ^ 1 

[12] [3670] = [123670] ~ 0 



coefficients 

products of matrices 

group index 

E 2 2c\ 

[65] J 8 = [65] ~ 0 


-E 2 2klcY 

[65] [2460] = -[5240] 

8 

E 2 2k 2 t 2 c ^ 1 

[65] [2470] = [652470] ~ 0 


-E^k^c^ 1 

[65] [2670] = [5270] 

5 

E 2 2c 2 k 2 Ci 1 

[65] [2463] = -[5243] 

10 

E 2 2c 2 t 2 c^ 1 

[65] [2473] = [652473] ~ 0 


E^clc^ 1 

[65] [2673] = [5273] 

3 

1 

[65] [3460] = -[5340] 

1 

E 2 2t 2 t 2 c ^ 1 

[65] [3470] = [653470] ~ 0 


—£/22c2f2Ci 1 

[65] [3670] = [5370] 

7 


There are 10 distinct groups of common terms in the previous eight tables. All the common terms cancel 
each other as shown by the following tables. 
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The 1st group: 


coefficients 

product of matrices 

2 fcici2ci 

[1672] 

—2fcici2&2hj" 1 

-[5340] = -[1672] 

because [1670] = [5243] [1672] = -[5043] 

-2tiC 1 2k 2 t2C[ 1 

-[5340] 

Ei2clc± l 

-[1267] = -[1672] 

—Ei2t 2 t 2 c^ 1 

[5340] 

—E 2 2k 2 t 2 c^ 1 

-[5340] 


The 2nd group: 


coefficients 

product of matrices 

— 2kiCi2k,2C^ 1 

[1740] 

2k\C\2c\ 

[2563] = -[1740] 

because [1670] [5243] [1470] = -[5263] 

-2tici2k 2 t 2 c^ 1 

[1740] 

— Ei2t 2 t 2 c^ 1 

[1470] = -[1740] 

E^c^c^ 1 

[5263] = [1740] 

E^k^c-^ 1 

[1470] 


The 3rd group: 


coefficients 

product of matrices 

2 c\2k2t2C 1 

-[1640] 

2tic\2ci 

[2573] = [1640] 

because [1670] = [5243] [1640] = -[5273] 

-2£ici2t 2 £2<h 1 

-[1640] 

—E^ktf^ 1 

[1460] = -[1640] 

E 2 2klcY 

[1460] 

E&clc^ 1 

[5273] - -[1640] 


( 107 ) 


(108) 


(109) 
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The 4th group of terms in the sum is 


coefficients 

product of matrices 

2k x c x 2c 2 k 2 E[ 1 

[1743] 

—2t\C\2c 2 t 2 c ^ 1 

-[2560] = -[1743] 

because [1670] = [5243] [1473] = [5260] 

-E^k^c^ 1 

[5260] = -[1743] 

—E 2 2c 2 t 2 Ci 1 

[1473] = -[1743] 


The 5th group: 


coefficients 

product of matrices 

2kici2c 2 t 2 c^ 1 

-[1643] 

—2kici2c 2 t 2 c± 1 

[2570] = -[1643] 

because [1670] = [5243] [1643] = [5270] 

—E 2 2c 2 k 2 c^ 1 

[1463] = -[1643] 

— E 2 2k 2 c 2 c 1 1 

[5270] = [1643] 


The 6th group: 


coefficients 

product of matrices 

— 2kiCi2c 2 t 2 Ci 1 

-[1230] 

2kici2c 2 t 2 c^ 1 

-[5647] = -[1230] 

because [1670] = [5243] ^ [5674] = [1203] 

2t 1 c 1 2c 2 k 2 c^ 1 

-[5746] = [1230] 

—2tici2k 2 c 2 c[ l 

-[1320] [1230] 


The 7th group: 


coefficients 

product of matrices 

—2fcic ;L 2fc 2 c 2 c^ 1 

[5370] 

2 £ici2c 2 t2<h 1 

-[1624] = -[5370] 

because [1670] = [5243] [5370] - [1246] 

^i2c 2 fc 2 c]f 1 

-[1246] = -[5370] 

—E 2 2c 2 t 2 c x 1 

[5370] 


31 


( 110 ) 


( 111 ) 


( 112 ) 


(113) 
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The 8th group: 


coefficients 

product of matrices 

-2k\C\2k 2 t 2 c{ x 

-[2540] 

2t\C\2t 2 t 2 c -^ 1 

-[2540] 

2t\C\2ci 

[1673] = [2540] 

because [1670] = [5243] [1673] = -[5240] 

—E\2k 2 t 2 c^ 1 

[5240] = -[2540] 

-E&clc^ 1 

[1673] 

-E 2 2k\c j -1 

-[5240] 


The 9th group: 


coefficients 

product of matrices 

-2fici2fc 2 c 2 cf 1 

-[5360] 

2t\ci2c 2 k 2 c^ 1 

[1724] = -[5360] 

because [1670] = [5243] 5=> [1274] = [5603] 

Ei2c 2 t 2 Ci 1 

-[1247] = [5360] 

—£^i2c 2 f 2 c 1 1 

[5360] 


The 10th group: 


coefficients 

product of matrices 

—£^i2c 2 f 2 c 1 1 
E\2c 2 t 2 c^ 1 
E 2 2k 2 c 2 c{ x 
E 2 2c 2 k 2 c x 1 

[1670] 

[5243] = [1670] 

[1670] 

-[5243] = -[1670] 


( 114 ) 


(115) 


(116) 
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